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Based on the complementarity relation between entanglement of a composite system and the 
purity of a subsystem, we propose a simple method to measure the amount of entanglement. The 
method can be applied to a bipartite system in a pure state of any arbitrary dimension. It requires 
only single qudit rotations and straightforward probability measurements performed on one of the 
subsystems, and can thus be easily implemented experimentally using linear optical devices. 
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One of the key issues in quantum information science 
is how to detect and measure entanglement. Perhaps the 
most straightforward way of measuring entanglement is 
full tomographic reconstruction of the quantum state in 
question. This technique is, however, highly inefficient 
and time-consuming. In recent years several less demand- 
ing methods of detecting and measuring entanglement 
have been proposed [1-7] and some demonstrated exper- 
imentally [8-10]. These methods typically require mea- 
surements performed on the entire composite system and 
involve coincidence measurements and/or complex con- 
trolled operations. Some methods require measurements 
upon two copies of the quantum state. 

On the other hand, it has been known that the amount 
of entanglement manifests itself in certain properties, 
such as the purity, of a subsystem [4,11,12]. Thus, the 
method of minimal and optimal tomography [13] on one 
of the subsystems can be used to measure entanglement 
[14] . Simpler ways of measuring the purity of a subsystem 
and consequently the entanglement exist [15,16]. These 
methods, however, are not without experimental difficul- 
ties. The method proposed by Ekert el al. [15] requires a 
controlled SWAP operation for which two identical copies 
of the state, whose purity is to be determined, need to 
be provided as target qudits. The method of Diirr [16] 
requires probability measurements with multiport beam 
splitters, in which probabilities need to be measured by 
varying phases associated with each port of the beam 
splitters over the entire range of the phase angles. 

In this paper we propose a simple method of measur- 
ing bipartite entanglement. The quantity to be measured 
is the concurrence [17] or its generalized version known 
as the / concurrence [11]. The system whose entangle- 
ment is to be measured needs to be in a pure state. The 
proposed method has several notable features. First, it 
requires measurements to be made only on one of the sub- 
systems. Second, it requires only single qudit rotations 
and straightforward projection measurements. Third, it 
can be applied to systems of any arbitrary dimension. 
We also describe how our method can be experimentally 
implemented using linear optical devices. 

The starting point of our analysis is the complementar- 
ity relation [11,12,18] between the / concurrence Cab of 
a composite quantum system AB and the purity Tr^(p^) 
of a subsystem A (or B), which reads, for the case when 



the composite system is in a pure state, 

Cis = 2 [1 - Tr^(p^)] = 2 [1 - TrsipD] (1) 

The amount of entanglement can thus be determined if 
the purity Tr(p^) of one of the subsystems is measured. 
As the purity can be written as 



Pij\ 



(2) 



where d is the dimension of the subsystem being mea- 
sured, d — dA (or d = ds) [From now on, we drop the 
subscript A (or B). It is to be understood that all quanti- 
ties refer to the subsystem being measured.], and 



means summation over all possible 



d(d-l) 



pairs of the 



indices i and j with i < j, the determination of the / 
concurrence requires knowledge on both diagonal and off- 
diagonal elements of the density matrix of the subsystem. 
The diagonal element pa is identified as the probability 
Pi of finding the subsystem in state \i) and can thus be 
determined directly from the projection measurement. 



Pi Pa ; 



i = 1,2, - • • ,d. 



(3) 



In order to determine the off-diagonal elements, we 
first need to perform a 90° rotation of the state of the 
subsystem about x axis in the clockwise direction in 
the three-dimensional Bloch sphere representation of the 
two-dimensional \k) — \l) space {k < I; fc, Z = 1, 2, • • • , d). 
The density matrix of the rotated system is given by 
p' = U[k,l)pU{k,l)\ where U{k,l) is a d x d rotation 
matrix with the elements given by 

1 i 

Ukk = Uii = — ^, Uki = Uik = 



Unn = 1 (n 7^ k,l), 



Unm — {n m; nm ^ kl, Ik). 

(4) 



The projection measurement performed on the rotated 
system yields the probabilities 

Pk = P'kk = \ [Pkk + Pii - i{pki - Pik)] , 



P! 
P' 



1 



Pii 

' Pnn ~ Pnn] 



2 [Pkk + Pll + i{pki - Plk)] , 



n ^ k,l. 



(5) 
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Since we know the diagonal elements pkk and pu from 
the straight projection measurement [Eq.(3)], we can de- 
termine the imaginary part of the off-diagonal element 
Pku Im(pfe;) = ^{pki - Pik), by measuring the probabil- 
ities Pf. and P/. This can be understood by noting that 
physically the 90° rotation about x axis projects the y 
component of the Bloch vector, which is equivalent to 
the imaginary part of the off-diagonal element, onto z 
axis. Performing the rotation U and the projection mea- 



surement on each of 



d(d-l) 



pairs of the states and 



1^), we can determine the imaginary part of all '^^'^^ 
off-diagonal elements pij. 

Having determined the imaginary part of the off- 
diagonal elements, we next need to determine the real 
part. That can be accomplished by a 90° rotation 
about y axis in the Bloch space representing the two 
states I A;) and The rotation brings the state to 
p" = V{k,l)pV{k,l)\ where 



Vkk 

Vnn = 1 (n ^ k,l), 



= 71' 



Vnm =0 [n^m; nm ^ kl, Ik). 



(6) 



The projection measurement on the rotated system yields 
now the probabilities 



Pk = Pkk = 2 IPkk + Pii + {Pkl + Plk)] 



Pi' = Pu = 2 \P^^ + ~ (^'^^ + ^'fc)] ' 

Pn ~ Pnn ~ Pnn] Tl ^ k,l. 



(7) 



The real part of the off-diagonal element pki, ILe{pki) = 
^{pki + Plk), can thus be determined from the measure- 
ment of the probabilities Pjl and P/'. The 90° rotation 
about y axis projects the x component of the Bloch vec- 
tor onto z axis, allowing the real part of pki to be de- 
termined. As before we need to perform the rotation V 
and the projection measurement on each of pairs 
of the states |fc) and |/) to determine the real part of all 
'^^'^2"^^ off-diagonal elements pij . 

Eq.(2) indicates that the determination of the purity 
requires the sum of the squares of all diagonal and off- 
diagonal elements, corresponding to the squares of the 
length of the Bloch vector projected onto z axis and onto 
xy plane, respectively. The sum of the squares of the 
diagonal elements can be obtained by squaring the prob- 
abilities Pi obtained from the straight projection mea- 
surement and summing them all. 



Information on the sum of the squares of the off-diagonal 



elements is contained in the quantity T' + T" where 

(k,l) i=l 



(k,l) i=l 



(9) 



Here Yli{k i) means summation over all possible ^^'^^ 
pairs of the states \k) and \l) with k < I. A straightfor- 
ward calculation yields 



T' + T" = {(f 



2d)T + 2J2 



+ 1. 



(10) 



Combining Eqs. (2), (8) and (9), we obtain for the purity 

Tr{p'^)=T' + T" -{(f -2d-l)T-l. (11) 

The / concurrence is then given according to Eq. (1) by 

Cab = V4 + 2(^2 -2d- 1)T - 2(T' -|- T"). (12) 

We now describe an experimental scheme that actu- 
ally performs the measurement. We adopt a multipath 
representation of a qudit. In a multipath system any ar- 
bitrary dx d unitary operator can be constructed using 
an arrangement of beam splitters, phase shifters and mir- 
rors [19]. In particular, the rotation operator U{k,l) of 
Eq. (4) can be realized by a 50/50 beam splitter situated 
at the intersection of the fcth and Zth paths as shown in 
Fig. 1 . The state transformation performed by the beam 
splitter is represented by the relation 



V2V 



^out 
^out 

'^out ~ ^inj 



(13) 



On the other hand, the rotation operator V{k, I) of Eq.(6) 
represented by the relation 



kout 

^out 



f^out 



nj^k 



(14) 



can be realized by a beam splitter at the intersection of 
the fcth and Ith paths and two phase shifters that shift 
phase by — ^ and |, respectively, situated at the input 
port and the output port, respectively of the Ith path, as 
shown in Fig.l. 

Our experiment then proceeds as follows. Wc assmne 
that a source emits a large number of pairs of path- 
entangled particles AB of dimension dA x ds = d x ds 
(see Fig.l). The / concurrence of the composite sys- 
tem AB is to be determined by performing local mea- 
surements on one of the subsystems, e.g., the subsystem 
A. The subsystem being measured has the dimension of 
dA = d. 
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FIG. 1: Rotation U(k,l)[V(k,l)] accomplished by a beam 
splitter (BS) placed at the intersection of the kth and Ith 
paths [and a — ^ phase shifter (— ) and a ^ phase shifter (+) 
at the input port and the output port, respectively, of the lih 
path]. 



(1) The projection measurement is performed directly on 
the subsystem A to measure the probability Pi for the 
subsystem A to be in the zth path. The measurement 
can be achieved with the arrangement shown in Fig.l 
with the beam splitter and the phase shifters taken out. 

Calculate T = X;ti 

(2) The projection measurement is performed with a 
beam splitter placed at the intersection of the fcth and 
Ith paths to determine the probability P/. The mea- 
surement can be achieved with the arrangement shown 
in Fig.l with the beam splitter kept installed but the 
phase shifters taken out. Repeat the projection mea- 
surement with the beam splitter shifted to the inter- 
section of all other possible pairs of paths. Calculate 

(3) The projection measurement is performed with a 
beam splitter and two phase shifters of — ^ and re- 
spectively, placed at the intersection of the fcth and llh 
paths to determine the probability P/'. The measure- 
ment can be achieved with the arrangement shown in 
Fig.l with the beam splitter and the phase shifters kept 
installed. Repeat the projection measurement with the 
beam splitter and the two phase shifters shifted to the 
intersection of all other possible pairs of paths. Calcu- 

later' = E(M)i:ti^^"'- 

(4) Determine the / concurrence using Eq.(12). 

The number of experimental setups required to com- 



each for steps (2) and (3), totaling (P — d + \. This 
number can considerably be reduced by placing beam 
splitters at every possible intersection simultaneously 
for steps (2) and (3), which increases the number of 
the real and imaginary parts, respectively, of the off- 
diagonal elements that can be determined from one ex- 
perimental setup to I for even d and for odd 
d. Each experimental setup for steps (2) and (3), 
respectively, then produces state p' and p" given by 
p' = U{kM, In)--- U{h,h)pUih,liy - - - U{kN, InY and 

■ vih,h)pv{hjiy ■■■vikMjNV, 



I for even d and N = ^ for odd d. This 



p" = VikN,lN) 

where N 

way the total number of experimental setups required 
reduces to 2(i — 1 for even d and 2d + 1 for odd d. A 
straightforward calculation shows that the / concurrence 
is then given by 



v /4 -I- 2{d - 3)r - 2(r^ + T") for even d 
y^A + 2{d-l)T-2{T' + T") for odd d 



(15) 



We illustrate our scheme for the simple case of a two- 

qubit system in state \iP)ab = cos 9\01) ab + sm9\10)AB ■ 

Let us choose the subsystem A to be the subsystem 

to be measured. Since p = pA — cos^6'|0)(0| -I- 

2/,i-,x/ii fcos^O \ ^. ,. ^ 
" 5> „ , the direct pro lection mea- 

sm^ 91' ^ 



•0|1)(1| 



surement should yield Pi = cos 9 and P2 ~ sin 9, and 
thus T = cos^ 9+8111^ 9. (Here the subscripts 1 and 2 refer 
to states |0) and |1), respectively.) It can be easily seen 

that p' = C/(l,2)p[/(l,2)t ^ ' ' ^ -ICOS29 



and 



p" = V{l,2)pV{l,2)^ 



Thus, 



i cos 29 1 
1 - cos 26*^ 
- cos 29 1 

the projection measurements of steps (2) and (3) should 
yield P{ = i, P^ = i and P\' = i, P^' = i, and there- 
fore T' = T" = i . The / concurrence for this case is the 
usual concurrence and is given according to Eq.(12) [or 
Eq. (15)] by Cab = \ sin 26*]. If the qubit is a photon po- 
larization qubit, the direct projection measurement cor- 
responds to a polarization measurement in the horizon- 
tal/vertical basis, while the projection measurement fol- 
lowing the rotation U and V, respectively, corresponds 
to a polarization measurement in the 45°/-45° basis and 
in the right-circular/left-circular basis. 

For the general two-qubit state |i/')a_b = a\00) ab + 
b\01)AB + c\10)ab + d\ll)AB, we obtain, through a 

straightforward calculation, Pi = [a'^ ' "^'^ 
12 p, _ 1 

_ 1 

Eq.(12) then immediately yields Cab 



\d\' 
Pi' 



2 + Im(ac* +bd*), P^ 
Re(ac* + bd*) and P^' 



plete the above experiment is 1 for step (1) and 



d{d+l) 



- Im(ac* + bd*), 

- Re(ac* +bd*). 
2\ad-bc\. 

In summary we propose a simple experimental scheme 
that allows one to determine the amount of entanglement 
of a bipartite composite system in a pure state. The 
scheme requires local probability measurements only on 
one of the subsystems. Furthermore, the scheme requires 
only single qudit rotations with linear optical devices and 
straightforward projection measurements. No complex 



4 



controlled operations are needed. The method can be 
applied to systems of any arbitrary dimension x cIb- 
The number of required experimental setups scales lin- 
early with the dimension of the subsystem being mea- 
sured. The scheme can be effective especially when one 
of the subsystems has a small dimension. If the compos- 



ite system is in a mixed state, Eq. (1) takes the form of 
inequality and our proposed method can only determine 
the upper bound of the I concurrence. 
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